Analysis of bulk and surface contributions in the neutron skin of nuclei 
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The neutron skin thickness of nuclei is a sensitive probe of the nuclear symmetry energy and has 
multiple implications for nuclear and astrophysical studies. However, precision measurements of this 
observable are difficult. The analysis of the experimental data may imply some assumptions about 
the bulk or surface nature of the formation of the neutron skin. Here we study the bulk or surface 
character of neutron skins of nuclei following from calculations with Gogny, Skyrme, and covariant 
nuclear mean-field interactions. These interactions are successful in describing nuclear charge radii 
and binding energies but predict different values for neutron skins. We perform the study by fitting 
two-parameter Fermi distributions to the calculated self-consistent neutron and proton densities. 
We note that the equivalent sharp radius is a more suitable reference quantity than the half-density 
radius parameter of the Fermi distributions to discern between the bulk and surface contributions 
in neutron skins. We present calculations for nuclei in the stability valley and for the isotopic chains 
of Sn and Pb. 
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Keywords: neutron skin thickness, halo-type nucleus, skin-type nucleus, antiprotonic atoms, halo factor, 
two-parameter Fermi distribution 
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I. INTRODUCTION 

The description of the sizes and shapes of atomic nu- 
clei is among the oldest problems in nuclear physics. The 
rms radius of the charge distribution in nuclei surely is 
the most prominent example of this kind of observable. 
Owing to the high degree of accuracy achieved by the 
elastic electron-nucleus and muon-nucleus scattering ex- 
periments, the nuclear charge radius is nowadays known 
with uncertainties that are for many nuclei smaller than 
1% [HQ. I n contrast, our knowledge about the neutron 
distribution and its rms radius in nuclei is not so precise. 
Actually, accurate determinations of the rms radius of 
the neutron density are still lacking. This implies that 
the so-called neutron skin thickness, generally defined as 
the neutron-proton rms radius difference in the atomic 
nucleus, 



Ai?„ p = 



,2x1/2 



2\1/ 



(1) 



is not precisely known either. 

The neutron skin thickness observable, Eq. (JTJ, is a 
very sensitive probe of the pressure difference that exists 
between neutrons and protons in the atomic nucleus. As 
such, AR np is intimately correlated with the density de- 
pendence of the nuclear symmetry energy and with the 
equation of state of pure neutron matter |3H12| . Owing to 
this fact, the accurate calibration of AR np is a problem of 
significant implications for studies that embrace diverse 
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facets of both nuclear phy sics and nuclear astrophysics 
(see for example Refs. [13l - l23l |). It should be mentioned 
that Eq. ([TJ is not the only useful prescription to charac- 
terize the different spatial extension of the neutron and 
proton densities in a nucleus. The neutron-proton radius 
difference has also been computed using Helm radii in- 
stead of rms radii for the discussion of nuclear halos in 
the literature [13, US]- Nevertheless, in the present work 
we use the conventional and more frequent definition ([T]) 
for ARnp, Eq. ©. 

Because neutrons are uncharged particles, the mea- 
surement of their spatial distribution in the nucleus is 
more difficult than for the positively charged protons 
(see, e.g., Ref. [26]). The experimental access to neutron 
densities and neutron rms radii usually involves strongly 
interacting hadronic probes, for example, in the case of 
experiments that perform proton-nucleus elastic scatter- 
ing d^-HU, cv-particle elastic scattering |32j ]. and tech- 
niques based on the inelastic scattering excitation of the 
giant dipole and spin-dipole resonances [33, 34] . The neu- 
tron skin thickness of nuclei has also been investigated 
through radiochemical and x-ray techniques in antipro- 
tonic atoms (35Uiol ] , taking advantage of the fact that the 
nuclear periphery is very sensitive to antiprotons in the 
normally electronic shell. Though the scope of our anal- 
ysis of neutron skins is of theoretical, in this article we 
refer to some extent to the experimental investigations in 
antiprotonic atoms to set the stage for our calculations. 

A few years ago, Trzcihska et al. [H, [3j| extracted the 
neutron skin thickness of a large set of nuclei in experi- 
ments with antiprotons conducted at the former LEAR 
facility of CERN. The measurements were made for 26 
stable isotopes distributed across the mass table, from 
the light and symmetric nucleus 40 C a to the heavy and 
asymmetric nucleus 238 U. Because of the fact that the 
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antiproton-nucleon interaction is very strong, antiprotons 
are able to interact with the atomic nucleus at distances 
where the nuclear density is much smaller than its central 
value. Slow enough antiprotons can form a hydrogen-like 
atom with the nucleus. When the antiproton annihilates 
with a nucleon producing pions that may miss the nu- 
cleus, it leaves a residue that is one neutron or proton 
fewer. From the analysis of these yields, information 
about the neutron distribution in the nucleus can be ob- 
tained (35l - l39| . A second experimental method measures 
antiprotonic x-rays from where the atomic level widths 
and shifts owing to the strong interaction are determined 
[38M40I ] . By combining the results obtained with these 
two experimental techniques, the neutron-proton rms dif- 
ference AR np can be deduced provided that the charge 



density of the nucleus is known [38H4 1 

The extraction of A.R np values from antiprotonic 
atoms assumes nucleon densities in the form of two- 
parameter Fermi (2pF) distributions. The procedure in- 
volves interpreting whether the difference between the 
peripheral neutron and proton densities arises from an 
increase of the mean location of the surface of the neu- 
tron density (i.e., from an increase of the bulk radius 
of neutrons) or, rather, from an increase of the sur- 
face diffuseness of the neutron density. This question 
is also instrumental in studies of properties of nuclei by 
parity- violating electron scattering [4J| . The radiochem- 
ical data in antiprotonic atoms were shown to be in favor 
of interpreting AR np as an increase of the neutron sur- 
face diffuseness [H, but with the caveat that some 
room existed within assigned errors for an intermediate 
situation 38[. Indeed, a comparison with the droplet 
model [42| suggests that the neutron skin sizes of the an- 
tiprotonic measurements can be described similarly well 
in the droplet model theory by a difference in the diffuse- 
ness as by a difference in the bulk radius of the neutron 
and proton densities. 

In the present work we theoretically investigate the 
bulk and surface components in the neutron skin thick- 
ness of nuclei by parametrizing self-consistently calcu- 
lated nucleon densities by 2pF distributions. The 2pF 
form is also common in experimental analyses. Our cal- 
culations are performed with some representative effec- 
tive nuclear forces. These are the finite-range Gogny D1S 
interaction and the zero-range Skyrme SLy4 force from 
the nonrelativistic framework, and the NL3 and FSUG- 
old parameter sets from the relativistic mean-field (RMF) 
framework. We discuss the predictions that these mean- 
field models make for the decomposition of the neutron 
skin thickness in bulk and surface contributions for the 
set of stable nuclei that were analyzed in the experiments 
in antiprotonic atoms. We also study the theoretical pre- 
dictions along the isotopic chains of Sn and Pb and the 
variation of the results as one moves from the proton to 
the neutron drip line in these isotopic chains. 

The structure of this article is the following. In the 
Sec. II we present a brief summary of the experimental 
methodology and results in antiprotonic atoms to high- 



light some interesting aspects for our study. In Sec. Ill 
we discuss the common definitions of nuclear radii and 
characterize the bulk and surface contributions in the 
neutron skin thickness of nuclei. In Sec. IV we ana- 
lyze the theoretical mean-field results in stable isotopes 
and the predictions for nuclei across the Sn and Pb iso- 
topic chains. We present the summary and conclusions 
in Sec. V. Some relations for 2pF functions are collected 
in the Appendix. 



SOME ASPECTS OF THE METHODS AND 
RESULTS IN ANTIPROTONIC ATOMS 



II. 



The radiochemical study of antiprotonic atoms [3E 
[39l . l43l - |45j consists of the analysis of the nuclei with a 
mass number one unit smaller than the target mass num- 
ber A t where the antiproton annihilation takes place. 
These products that have one less neutron or proton are 
short lived and their decay is followed by emission of 7 
rays. Standard nuclear spectroscopy methods allow one 
to determine the absolute number of these residual nuclei 
with mass number At — 1. The probability of producing 
a cold product with mass number A t — 1 is calculated 
by using the antiproton-nucleus optical potential fitted 
to reproduce the x-ray experimental data (atomic level 
shifts and level widths) [H, H(| Efl • The probability dis- 
tribution for obtaining a nucleus with mass number At — 1 
has a maximum located about 2-3 fm outside the half- 
density radius R1/2 of the target nucleus. 

To compare the experimental data for any target it is 
convenient to introduce the peripheral halo factor 471 ] : 
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The first term in brackets on the right-hand side of this 
equation, gives the ratio of the p annihilations on pe- 
ripheral neutrons to the p annihilations on peripheral 
protons, normalized with the Z t /N t value of the target 
nucleus. The quantities Sp n and Sp p are the p-n and 
p-p scattering amplitudes, respectively, and the factor 
F = Im(<5p n )/Im((5pp) gives the ratio of the p annihila- 
tion probabilities on a neutron and on a proton. Conse- 
quently, halo factor ([2]) essentially measures the change 
of the ratio of the neutron-to-proton concentration in the 
peripheral region with respect to the bulk value repre- 
sented by the N t /Z t ratio in the target nucleus. A halo 
factor larger (smaller) than 1 means an increase of the 
relative neutron (proton) concentration in the nuclear 
periphery. Assuming the ratio F is known, the measure- 
ment of the N t — 1 and Z t — 1 product yields allows one to 
obtain the experimental value of the halo factor. It was 
concluded that the comparison of the results from the an- 
tiprotonic x-ray method with the radiochemical method 
favors a value of F m 1 [38| . 

The number of nuclei withone fewer neutron or proton 
is proportional to the quantities lva(8 pn )p n or lm(5 pp )p p , 
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respectively, integrated over the region where the annihi- 
lation process takes place. Thus, approximately, one can 
estimate theoretically the halo factor as [33, 0] 
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Pn{r) Z t 
P P (r) N t 



(3) 



This factor reaches the experimental value f^f Q at a dis- 
tance r « R1/2 +2.5 fm beyond the half-density radius 
of the charge distribution, where the antiproton annihi- 
lation takes place. One can use this fact to reproduce 
the neutron density distribution from some experimental 
observables. 

The charge distribution in many stable nuclei is known 
very precisely. Usually the experimental charge density is 
given in some analytical form that is or may be converted 
to a 2pF distribution [H, H3| • The simple 2pF formula 



p(r) = 



po 



1 + cxp [(r - C)/a] 



(4) 



has only two free parameters with clear physical mean- 
ing: on the one hand, a describes the diffuseness of the 
surface of the density profile: on the other hand, the half 
density or central radius C describes the mean location 
of this surface (i.e., C is indicative of the extension of the 
bulk part of the density distribution). The other, more 
sophisticated, parametrizations of the density distribu- 
tions [4!| describe better the central part of the nu- 
cleus or modify the surface part with higher-order terms. 
We do not use them because the interpretation of the 
multiple parameters is harder and less direct. 

To study the differences at the nuclear periphery be- 
tween the neutron and proton densities in the antipro- 
tonic atoms experiments, the authors of Ref. [38[ used 
2pF functions. They used the notation "neutron skin- 
type" distribution and "neutron halo-type" distribution 
to describe the two extreme cases of 2pF shapes having 
either C„ > C p and a n — a p , or C n = C p and a n > a p , 
respectively. As we mentioned in the Introduction, the 
radiochemical data gave support to understanding the 
difference between the neutron and proton rms radii in 
the antiprotonic atoms from an increase of a n rather than 
from an increase of C n , compared to the a p and C p values. 
Thus, the neutron halo-type distribution (C„ — C p = 0) 
was assumed [38| . (Note that here "halo type" is a useful 
notation, but the nuclei in the antiprotonic experiments 
are stable isotopes and it is not meant that they have ha- 
los like very light or exotic nuclei, such as 11 Li 0, 
It was noted in the same work [38| that some room was 
left, nevertheless, within the error bars for intermediate 
cases with C n > C p and a n > a p . 

The 2pF shape can be applied for the description of 
charge, proton, or neutron densities. Indeed, if the charge 
density is known in the 2pF form, the corresponding 
point proton density can be easily found and it also takes 
a 2pF form, with the parameters obtained through the 
deconvolution procedure [5lL HH . The more relevant ex- 
pressions are given in the Appendix. From these proton 



DIS 

SLy4 

C„-C p = -0.20 - 

C„-C p = -0.10 

C„-C p = 0.00 - 

C„-C p = 0.10 

C„-C p = 0.20 - 


,///, ' ' 

t/Af s 

"/ 


^'^ 





6 

r(fm) 



10 



12 



FIG. 1: (Color online) The halo factor in 208 Pb as a func- 
tion of the distance from the center of the nucleus, calculated 
usin g th e 2pF formulas from the experimental charge den- 
sity [53] assuming a neutron skin thickness AR np — 0.16 fm 
and different values of C'„ — C p (thin lines). The results of the 
theoretical predictions with the Gogny DIS and Skyrme SLy4 
nuclear forces are also shown. The values of the halo factor 
deduced from experiment [38|, are marked by crosses. The 
value of the proton half-density radius C p is indicated by an 
arrow. 



density profiles one can try to deduce the neutron density 
profiles with some additional information. For instance, 
if one knows the C p and a p parameters of the proton den- 
sity and the value of the neutron skin thickness AR np 
defined in Eq. ([IJ (e.g., from independent experimental 
measurements), a relationship between the half-density 
radius C n and the diffuseness a n of the neutron 2pF dis- 
tribution can be found (see the Appendix). Such a pro- 
cedure creates a family of neutron density profiles that 
depend on one free parameter, which can be chosen to 
be either C n — C p or a n — a p . Now it is possible to check 
which set of values of C n and a n reproduces the experi- 
mental halo factor ([2]) . This allows one to determine the 
2pF neutron distribution. Next we apply this idea to the 
208 Pb nucleus as an illustrative test case. 

In Fig. Q] we display the halo factor, Eq. ([3]), for the 
208 Pb nucleus computed with 2pF density distributions. 
The 2pF neutron densities are obtained from the exper- 
imental charge density taken from Ref. [53j, which we 
convert to a 2pF proton density using the formulas given 
in the Appendix. We consider a family of 2pF neutron 
densities that differ one from the other in the half-density 
radius C n and diffuseness a n , but they all have the same 
value for the neutron skin thickness: AR np — 0.16 fm 
(which corresponds to the value extracted from the re- 
sults of the x-ray method in the 208 Pb antiprotonic atom 
[40|). We allow the difference C n — C p to vary in the 
range from —0.20 to 0.20 fm. The value of the diffuse- 
ness a n of the 2pF neutron density is then obtained using 
Eq. (|A10p from the Appendix. 

The values of the 2pF parameters in our present cal- 
culation are given in Table HI Agreement with the results 
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TABLE I: The parameters C and a of the 2pF distributions 
for 208 Pb obtained from the experimental charge density [53| 
as described in the text and from mean-field calculations. All 
values are given in femtometers. 





C g 


a q 


C n Cp 


0"TI 


proton 


6.704 


0.438 






neutron 


6.504 


0.659 


-0.20 


0.221 




6.604 


0.614 


-0.10 


0.176 




6.704 


0.565 


0.00 


0.127 




6.804 


0.511 


0.10 


0.073 




6.904 


0.449 


0.20 


0.011 


fit to D1S density profile 


proton 


6.645 


0.467 






neutron 


6.686 


0.548 


0.041 


0.081 


fit to SLy4 density profile 


proton 


6.683 


0.470 






neutron 


6.755 


0.555 


0.072 


0.085 



from previous experiments [38l |40j (indicated by crosses 
in Fig. [T]) is found for the parameters of the 2pF distri- 
butions in the range from the values C n — C p = fm and 
a n — a p = 0.13 fm to the values C n — C p — 0.10 fm and 
o„ — a p — 0.07 fm. In the first limit, the neutron skin 
thickness AR np is basically due to an enhancement of the 
diffuseness of the neutron density, whereas the second 
limit corresponds to a mixed configuration where both 
the neutron diffuseness and the neutron half-density ra- 
dius are larger than the values of the proton density. We 
have to stress that even small differences between C n and 
Cp may have a meaningful influence on the value of the 
neutron skin thickness. The discussed results indicate 
an intermediate character of the neutron skin thickness 
in 208 Pb, but with a preference for the density pattern 
where C n ~ C p and a n > a p . A more precise conclusion 
is difficult because of the large experimental error bars. 

For comparison, in Fig. Q] we also plot the theoreti- 
cal halo factor ^ in 208 Pb calculated directly from the 
neutron and proton densities of the Gogny D1S |54| and 
Skyrme SLy4 [55| effective nuclear forces. At the dis- 
tances relevant for antiproton annihilation, the halo fac- 
tor obtained with these theoretical mean-field densities 
agrees considerably well with the experimental data. The 
numerical values of the half-density radius C and of the 
diffuseness parameter a of the D1S and SLy4 equivalent 
2pF neutron and proton densities for 208 Pb are reported 
in Table [I] For either force, one observes that the half- 
density radii C n and C p of the 2pF distributions are dif- 
ferent and that the values of a n and a p are also different. 
Thus, the results of these models correspond to some 
intermediate situation between the "neutron halo-type" 
2pF distribution (where C n — C p and a n > a p ) and the 
"neutron skin-type" 2pF distribution (where C n > C p 
and a n — a p ) discussed in Ref. (38]. In any case, the 
two forces, D1S and SLy4, show some preference, as 
also do the calculations considered earlier, for the "neu- 
tron halo-type" distribution in 208 Pb, especially in the 
case of the D1S force, where C n — C p = 0.04 fm and 



a„ — a p — 0.08 fm. 

From the discussions of this section it is easily seen that 
the characterization of the "bulk" or "surface" forma- 
tion of the neutron skin thickness in nuclei is a nontrivial 
problem. In the following, we wish to analyze this ques- 
tion according to the predictions derived from mean-field 
models of nuclear structure that are tested to be suc- 
cessful for charge radii, binding energies, and a wealth 
of phenomena in nuclei. First, we need to establish a 
prescription to separate bulk and surface contributions 
in neutron skins calculated with mean-field nuclear den- 
sities. 



III. DISCERNING BULK FROM SURFACE IN 
NEUTRON SKINS 

Nuclear radii can be defined in several independent 
ways. The most popular formulas and their relations 
are discussed in the book by Hasse and Myers [49[ and 
we recall and analyze them in this section. One of the 
simplest ways to describe the size of nuclei is to define 
a central radius C in terms of the integral of the nuclear 
density profile p{r) (for either neutrons or protons) as 

c = W)[* r)dr - (5) 

Another option is the half- density radius R\/2 which is 
defined from the local condition 

p(Ri/2) = \ P {0). (6) 

For density profiles that have a symmetric surface, such 
as the 2pF distribution defined in Eq. (|4]), the central 
radius C and the half-density radius R\/2 coincide. 

The equivalent sharp radius R is the radius of a uniform 
sharp distribution that has a constant density equal to 
the bulk value of the actual density and that contains the 
same number of nucleons as the considered nucleus: 

4 f°° 

-irR 3 p{b\Ak) = 4tt / p{r)r 2 dr. (7) 
3 Jo 

For its importance in experimental techniques, the equiv- 
alent rras radius Q is commonly used. It describes a uni- 
form sharp distribution with the same rms radius as the 
given density profile: 

| Q 2 = (r 2 ) • (8) 

Because the neutron skin thickness ([T]) is defined through 
rms radii, it can be expressed easily with Q: 

ARn P = J^(Q n -Q p ). (9) 

In uniform, sharp-edge nuclear distributions all of the 
above-mentioned definitions coincide, but in realistic lep- 
todermous density profiles they give distinct values. The 
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FIG. 2: (Color online) Comparison of sharp surface density 
profiles that have the central (C), equivalent sharp (R), and 
equivalent rms (Q) radii with the self-consistent and 2pF pro- 
files corresponding to the neutron density of (a) 100 Sn and (b) 
132 Sn obtained in the RMF theory (the NL3 interaction [5fj] 
has been used). 



relation between the C, Q, and R radii can be expressed 
by expansion formulas in powers of b/R, where b is the 
so-called surface width of the density profile. The latter 
quantity is defined by 



b 2 = - 



P(0)J ^ C) dr 



dr . 



To leading order, one has the relationships |4£ 
C~R 

and 



R 2 ) 



(10) 



(11) 



(12) 



These expansions are useful as long as b/R is small. This 
condition is well fulfilled in many nuclei because b ~ 1 fm 
and R ~ tqA 1 ^; therefore, b 2 / R 2 is typically no larger 
than A -2 / 3 . Moreover, to consider no further corrective 
terms in the above relations of the C and Q radii with R 
may be quite accurate for specific shapes. For example, 
in the case of 2pF distributions, the first nonvanishing 
corrections to the terms within brackets in Eqs. (fTTj) and 
(IT21 are of order (b/R) 6 and (b/R) 4 , respectively. 

The multiple definitions of nuclear radii may some- 
times cause misleading conclusions, especially if nuclear 



properties sensitive to the nuclear matter distribution in 
nuclei are concerned or two different models are com- 
pared. Therefore one has to be careful about the suitable 
choice of the radius definition. To compare the forego- 
ing definitions of radii for heavy nuclei, in Fig. [2] we plot 
the neutron densities obtained in a self-consistent mean- 
field calculation and the fitted 2pF distributions for the 
neutron-deficient nucleus 100 Sn and the neutron-rich nu- 
cleus 132 Sn. These profiles are compared with sharp-edge 
density distributions having radii C, R, and Q, calculated 
from the above expressions with the 2pF function. The 
central densities of the sharp surface spheres are fixed so 
to fulfill the particle number normalization. 

Figure [5] illustrates the fact that the central radius C 
does not allow the bulk density to be reproduced. A 
sharp sphere of radius C overestimates the self-consistent 
density in the whole nuclear interior. The equivalent rms 
radius Q also fails, because it clearly underestimates the 
original density in the bulk. Only the equivalent sharp 
radius R is able to properly reproduce the bulk part of 
the self-consistent and 2pF density profiles of the nucleus. 
As discussed in Ref . (4{| , the equivalent sharp radius R is 
the quantity of basic geometric importance of the three 
radii C, R, and Q. A sharp distribution of radius R 
has the same volume integral as the actual density of 
the finite nucleus and differs from it only in the surface 
region. Therefore, the radius R appears to be the suitable 
quantity to be used to measure the size of the bulk part 
of the nucleus. 

On account of Eq. © for the neutron skin thickness 
and relationship (fT2")) between Q and R, one obtains the 
expression 



AR np = \l- 



(R n - R p) + \\jr-jr 

Z \ K n rip 



(13) 



up to terms of order 0(b 4 /R 3 ). Thus, one can make a 
meaningful distinction between a bulk contribution and 
a surface (diffuseness) contribution to the neutron skin 
thickness of nuclei as follows: 



bulk 



with 



A psurl 



(R n — R p ) 



independent of surface properties and 

bl 



AR 



surf 




(14) 



(15) 



(16) 



The nucleus may develop a neutron skin by separation 
of the bulk radii R of neutrons and protons or by mod- 
ification of the diffuseness b of the neutron and proton 
surfaces. In the general case a combination of both ef- 
fects can be found. To which degree the different pat- 
terns arise in mean-field calculations of finite nuclei is a 
question we address in the following sections. 
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Experimental [50| and theoretical mean-field density 
distributions present oscillations in their inner bulk re- 
gion, which implies that some suitable average is needed 
to determine the bulk density value p(bulk) of Eq. 
to compute the equivalent sharp radius R. The difficulty 
may be easily solved by fitting a Fermi function to the 
original density, as we have illustrated in Fig. [5J In this 
case, the bulk density value ,that is, po/[l +cxp (— C/a)], 
is, to excellent accuracy the po parameter of the Fermi 
function. An effect not described by functions like the 
2pF distribution is the occurrence of a nonsymmctric 
surface shape, which is possible in real nuclear densi- 
ties. However, one plausibly expects that neutron skins 
of nuclei are dominated by the difference existing between 
neutrons and protons in the location of the surface (i.e., 
R) and/or in the spatial extent of this surface (i.e., 6); 
the difference in the degree of asymmetry between the 
shapes of the neutron and proton surfaces is expected to 
be a less important, higher-order correction. 

It may be practical to rewrite expressions (| 15|) and (|16|) 
for the bulk and surface contributions to the neutron skin 
thickness directly in terms of the parameters of the 2pF 
function of Eq. Q. Using the fact that the diffuseness 
parameter a in a 2pF function is related to the surface 
width b by the formula 



(17) 



and inverting relation (fTTj) between C and R, one easily 
finds to the given order of approximation that Eqs. (|15j) 
and p^l) become, respectively, 



Ai& lk = \l- 







4\ 






c p )_ 



(18) 



and 



■surf 




AR 



We note from these results that 



(19) 



AR 



bulk 

np 



^(C n ~C p ) + ^AR SUI 
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■up 



(20) 



that is, the quantities C n — C p (easily obtained from 
the 2pF distributions themselves) and R n — R p [obtained 
through Eq. ([?])] differ by surface diffuseness terms and 
should not be mixed. We have seen in Fig. [2] that a 
sharp sphere having a radius C significantly distorts the 
appearance of the actual nuclear density by overshooting 
it in the whole bulk region. It is thus preferable to use 



3/5(i?„ — R p ) rather than 



■ C„ 



as a measure 



of the bulk contribution to the neutron skin thickness. 



IV. DECOMPOSITION OF NEUTRON SKIN 
THICKNESS IN SELECTED NUCLEI - 
MEAN-FIELD RESULTS 

A. Nuclei of the experiments in antiprotonic atoms 

To get information about the "bulk" or "surface" 
character of the thickness of neutron skins in theoret- 
ical mean-field calculations, we parametrize the self- 
consistently calculated proton and neutron densities with 
2pF distributions. This procedure can be applied very 
suitably for many heavy nuclei and gives a clear distinc- 
tion between bulk and surface properties of nuclei. How- 
ever there is no universal method to do this parametriza- 
tion. A popular prescription is to use a x 2 minimization 
of the differences between the density to be reproduced 
and the 2pF profile, or of the differences between their 
logarithms. These methods may somewhat depend on 
conditions given during minimization (number of mesh 
points, limits, etc.). We prefer to extract the parameters 
of the 2pF profiles by imposing that they reproduce the 
same quadratic (r 2 ) and quartic (r 4 ) moments of the self- 
consistent mean-field densities. These two conditions, 
together with the normalization to the proton and neu- 
tron numbers, allow us to determine in a unique way the 
equivalent 2pF densities. This method can be applied to 
any density distribution. Its focus is on a good repro- 
duction of the surface region of the original density be- 
cause the local distributions of the quantities r 2 p(r) and 
r 4 p(r) are peaked at the periphery of the nucleus. As 
an example of the results of the present determination 
of the 2pF profiles, in Fig. [3] we display in logarithmic 
scale the self-consistent neutron and proton densities for 
208 Pb together with their equivalent 2pF distributions 
calculated with the NL3 interaction. It can be seen that 
there is overall good agreement in the central and surface 
regions of the nucleus. In particular, the 2pF densities 
reproduce the mean-field densities well at the distances 
that are relevant for antiproton annihilation. 

In this section, we would like to get some insight about 
the "bulk" or "surface" character of the neutron skin pre- 
dicted by the mean-field approach in the nuclei for which 
the neutron skin thickness values were extracted in Refs. 
[HI, [3^1 from the measurements in antiprotonic atoms. 
These nuclei range from 40 Ca to 238 U and all lie along 
the valley of stability. We carry out the calculations with 
the aforementioned nonrelativistic interactions D1S 1 54] 
and SLy4 [55[ plus the relativistic interactions NL3 [56 j 
and FSUGold [57| , as representative examples of success- 
ful nuclear mean-field models. We impose spherical sym- 
metry in all nuclei described in this article. The effect of 
deformation on the neutron skin thickness was discussed 
elsewhere [58[ . The two nonrelativistic forces have a soft 
symmetry energy (l5l . [l8j . On the contrary, the covariant 
NL3 parameter set has a stiff symmetry energy, which is 
usual in the relativistic models [l5|, [l8[ ■ Note that the 
notation soft or stiff refers to whether the symmetry en- 
ergy of the model increases slowly or rapidly as a func- 
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r(fm) 

FIG. 3: (Color online) The proton and neutron densities of 
208 Pb calculated with the NL3 parameter set as a function 
of the distance from the center of this nucleus in logarithmic 
scale. The 2pF densities fitted to the NL3 results are also 
shown. 



tion of the nuclear density. The covariant parameter set 
FSUGold was devised to have a softer density dependence 
of the symmetry energy (57l | than the typical relativistic 
models. Thus, FSUGold is, in this aspect, closer to the 
nonrelativistic models than NL3. 

In Fig. 2] we display the experimental data with error 
bars determined from the antiprotonic atoms. There is 
a relatively clear correlation between the experimental 
value of the neutron skin thickness of these 26 stable nu- 
clei and the overall relative neutron excess I = (N—Z)/A 
of the nucleus. This trend has been fitted by the linear 
relationship AR np = (0.90 ± 0.15)/ + (-0.03 ± 0.02) fm 
with a x 2 factor of 0.5 38, 39]. In the same figure we plot 
the theoretical AR np value calculated according to Eq. 
([1]) using the mean-field densities of the indicated effec- 
tive nuclear models for the 23 even-even nuclei that exist 
in the experimental data set. It is obvious that the the- 
oretical models make largely different predictions for the 
neutron skin thickness. This is especially visible at large 
values of /. It is seen that the models that have a softer 
symmetry energy give smaller AR np values, whereas the 
models with a stiffer symmetry energy give larger AR np 
values [U EH- Note that as soon as / ^ 0, even when 
it is small, discrepancies arise among the AR np values of 
the models. In contrast, all of the considered interactions 
make an almost identical prediction of AR np rj —0.05 fm 
for the 40 Ca nucleus. 

One observes in Fig. @] that, similarly to the situa- 
tion in the experimental data set, the theoretical neu- 
tron skin values computed with each nuclear interaction 
show an average linear behavior as a function of the 
neutron excess /. Actually, the linear correlation fac- 
tor of AR np with / in the present models is considerably 
high (between 0.95 and 0.97). The fit of the neutron 
skin values calculated with the Skyrme SLy4 force yields 
AR np = 1.01/ - 0.035 fm, whereas AR np = 0.88/ - 0.04 



0.3 1 — | 1 1 1 1 1 1 r 




J i Li i_l i I i i i i L_ 

0.1 0.2 
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FIG. 4: (Color online) The results of the covariant NL3 and 
FSUGold parameter sets and of the nonrelativistic Skyrme 
SLy4 and Gogny D1S forces compared with the experimen- 
tal neutron skin values AR np deduced from antiprotonic 
atoms (solid squares with errorbars) and their linear average 
A_R^= (0.90 ± 0.15)/ + (-0.03 ± 0.02) fm (shaded region) 



fm is obtained in the case of the Gogny D1S interac- 
tion. We find a linear fit of AR np = 1.20/ - 0.03 fm 
using the values of the neutron skin thickness computed 
with the relativistic FSUGold model. If we consider the 
relativistic NL3 parameter set, which has a stiffer sym- 
metry energy than the other models, we find the linear 
fit AR np — 1.50/ — 0.03 fm. Compared with the slope 
of 0.90 for the fit of the experimental data [lH, the 
slopes in forces like D1S and SLy4 is much closer to it, 
whereas the agreement deteriorates as the model has a 
stiffer symmetry energy. Thus, one can conclude that the 
comparison of the slopes of AR np with / between theory 
and the antiprotonic measurements for nuclei across the 
mass table favors the models that have a soft symmetry 
energy. 

In Fig. [5] we display the bulk and surface contribu- 
tions to the theoretical neutron skin thickness that are 
computed by applying Eqs. (fTS")) and (flB)) , respectively, 
against the neutron excess /. These values are obtained 
from the 2pF distributions associated with the mean-field 
neutron and proton densities calculated with the D1S 
force as well as with the Skyrme SLy4 force and the two 
RMF parameter sets FSUGold and NL3. The numeri- 
cal calculations show that the value of A/?^ lk + AR s ^ p f 
obtained through Eqs. (fT5| and (fTI)|) can be slightly less 
accurate in some nuclei [compared to the exact value of 
Eq. (JTJ>] than the result obtained through Eqs. ([15]) and 
(|16[) . The small differences, whenever they arise, are due 
almost entirely to the replacement of a q /R q by a? q jC q 
(q = n,p) in Eqs. flTSJ) and ([19]), Of course, Eqs. $TSj) and 
(|19[) are quite practical because they do not require the 
additional calculation of the R q values and can be ap- 
plied straightforwardly using the C q and a q parameters 
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FIG. 5: (Color online) The neutron skin thickness AR np as well as its bulk Ai?J^ lk and surface Ai?^p rf parts [given by Eqs. 
(|15[l and (|16|) . respectively] for the nuclei presented in Fig. [4] calculated in various models. 



of the 2pF distributions themselves. 

The four panels presented in Fig. [5j as a consequence 
of the differences between models, predict a slightly dif- 
ferent splitting of the neutron skin thickness into their 
bulk and surface contributions for each considered nu- 
cleus. Therefore, the values of the bulk and surface con- 
tributions to AR np are to some extent model dependent. 
The total neutron skin thickness roughly follows an in- 
creasing tendency with the relative neutron excess /, as 
discussed previously. This tendency is also seen in the 
bulk contribution A-R^" , especially in the nonrelativis- 
tic interactions. However the surface contribution does 
not clearly follow the same trend and shows a less definite 
behavior as a function of /. 

From Fig. [5] we see that for relatively neutron-rich nu- 
clei (/ > 0.15), the surface part generally contributes 
50% or more to the total neutron skin thickness. If the 
nuclear model has a stiff symmetry energy and / is large, 
the bulk part of AR np may become larger than the sur- 
face part, as can be seen in the NL3 panel of Fig. [5] 
at / > 0.2. More symmetric nuclei (/ < 0.15) do not 
present a definite tendency. The theoretical calculations 
show that in these nuclei the sharp radius is larger for 
protons than for neutrons (R p > R n ) while the surface 
width is larger for neutrons than for protons (b n > b p ). 
Consequently, the bulk contribution to the neutron skin 
thickness becomes negative, as it can be seen in Fig. [SJ 
and the relatively small value of the neutron skin thick- 
ness is basically because of a strong cancellation between 
the bulk and surface parts. In the lightest nuclei, both 
contributions are negative and produce a "proton skin" 
rather than a neutron skin. 



B. Medium and heavy mass isotopic chains 

The set of nuclei chosen in Sec. IV A have given us 
some insight into the bulk and surface contributions to 
the neutron skin thickness in stable isotopes. To investi- 
gate in a systematic way how the bulk and surface con- 
tributions evolve with the neutron number we study the 
chains of even-even Sn and Pb isotopes from the proton 
to the neutron drip line. 

First, Fig. [6] displays the neutron skin thickness AR np 
along the Sn and Pb isotopic chains computed in the 
mean-field approximation with the Gogny D1S and 
Skyrme SLy4 forces and using the relativistic mean-field 
models NL3 and FSUGold. For small and moderate val- 
ues of the relative neutron excess (/ < 0.2), the neutron 
skin thickness grows almost linearly with / in each iso- 
topic chain, as in the case of the stable nuclei analyzed 
in the previous section. However, AR np shows a rather 
pronounced kink at / ss 0.20 — 0.25. Beyond this value, 
it increases again almost linearly as a function of the rel- 
ative neutron excess, but with a larger slope. Finally, 
a new departure from the linear behavior of AR np as 
a function of / can be observed when the isotopes are 
on the edge of the neutron drip line. The kinks and 
changes of slope in the neutron skin thickness as a func- 
tion of the relative neutron excess are clearly connected 
with the doubly magic nuclei 132 Sn (7 a 0.24) and 208 Pb 
(/ a 0.21) in the stable nuclei region and with another 
two doubly magic nuclei, namely 176 Sn (/ ss 0.43) and 
266 Pb (/ « 0.38), near the neutron drip line. Therefore, 
it is obvious that these changes of slope are produced by 
quantal effects, which modify the linear trend of AR np 
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FIG. 6: (Color online) The neutron skin AR np for (a) Sn and 
(b) Pb isotopes calculated with several mean-field models. 



as a function of I in a non-negligible way. The average 
slope of AR np with / for various forces is clearly differ- 
ent. As has been shown in Refs.fTTI. fl2j. for each nuclear 
model, AR np is strongly correlated with the slope of the 
symmetry energy with respect to the density computed 
at saturation. 

As in Sec. IV A, we fit the mean-field proton and neu- 
tron densities by 2pF distributions to investigate the bulk 
and surface contributions to the neutron skin. To ana- 
lyze the changes that occur in AR np first we look at the 
parameters that characterize the 2pF distributions. In 
Fig. [7] we display the central density po, the half-density 
radius C and the diffuseness a, which are obtained from 
the D1S densities along the Sn and Pb isotopic chains. 
The 2pF parameters show an overall smooth behavior as 
a function of the relative neutron excess / but with local 
modulations near the shell closures. The central density 
Po grows (for neutrons) or declines (for protons) almost 
linearly as a function of / in both isotopic chains. It is 
just a simple consequence of the increasing asymmetry 
along the isotopic chains. The central radii C n and C p 
both show a global increasing tendency with the neutron 
number. One can notice that C n grows faster with / 
than C p . However, the a n and a p diffuseness parameters 
behave in a completely different way. Whereas a n grows 
on average with increasing neutron excess, a p remains 
roughly constant with / in the two isotopic chains. On 
top of these general trends we see that the 2pF param- 
eters associated with the neutron densities show kinks 
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FIG. 7: (Color online) Proton and neutron 2pF parameters 
(a) Po for Sn, (b) p for Pb, (c) C for Sn, (d) C for Pb, (e) 
a for Sn, and (f) a Pb isotopes. 2pF parameters are fitted to 
the Gogny D1S density distributions. 



around the shell closures at N = 82 in Sn and N = 126 
in Pb (that is, in the region / ~ 0.20 — 0.25), as well 
as changes of slope near the neutron drip lines. One ob- 
serves that some signature of the neutron shell effect also 
appears in the 2pF parameters of the proton densities 
around magic neutron numbers. 

Figure [S] shows the evolution with / of the 2pF param- 
eters corresponding to the quantal densities computed 
with the SLy4 Skyrme force and with the relativistic 
mean-field NL3 parameter set. In the two models the 
global trends are similar to the ones discussed before for 
2pF parameters of D1S distribution. However, they show 
some differences that can be attributed in part to the 
different properties of the isovector channel of the inter- 
action. 

In Figs.^a) and[SJb) we display the bulk and the sur- 
face contributions to the nuclear skin thickness calculated 
with the D1S force along the Sn and Pb isotopic chains. 
These contributions show for Sn isotopes two well-defined 
regions as a function of /. One of them covers the neutron 
major shell between 100 Sn and 132 Sn, i.e. < I < 0.25, 
and the other region corresponds to the next major shell 
between 132 Sn and 176 Sn in the range 0.25 < I < 0.43. 
For nearly symmetric Sn isotopes close to 100 Sn, which 
are neutron deficient, C p is larger than C n (see Fig. [7]) 
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FIG. 8: (Color online) The same as in Fig. [Jj but for the 
Skyrme SLy4 force and relativistic NL3 parameter set. 



and, consequently, the bulk part of the neutron skin is 
negative or at most it takes very small positive values. 
For these values of /, the surface contribution is posi- 
tive and relatively small, reducing the opposite effect of 
"proton skin" due to negative values of the bulk part. 
Looking at more asymmetric isotopes, one can see that, 
in the magic nuclei 132 Sn and 176 Sn, the bulk and the 
surface contributions are roughly equal. This implies a 
rather compact neutron density distribution with a rela- 
tively stiff surface as a consequence of the kinks exhibited 
by the neutron diffuseness parameter a n around the neu- 
tron magic numbers, which can be seen in Fig. [7J In 
the regions between magic numbers, the surface contri- 
bution to the neutron skin thickness is larger than the 
bulk contribution. The splitting between the surface and 
bulk contributions reaches its maximum value roughly at 
midshell. This behavior of the bulk and surface contri- 
butions to the neutron skin thickness points out that Sn 
isotopes in the middle of major shells develop a larger 
surface region in the density distribution than the magic 
ones. In other words, these isotopes with neutron num- 
ber in between magic values are more of "halo" type than 
the limiting magic nuclei which show a mixed character 
between "halo" and "neutron skin." 

The bulk and the surface contributions to the neutron 
skin of the Pb isotopes calculated with the D1S force 
show a similar behavior to the case of the Sn isotopes 



analyzed before. However, for this heavy isotopic chain, 
the bulk part gives a more important contribution to the 
total neutron skin for nuclei in between the magic 208 Pb 
and 266 Pb nuclei. 

The discussed differences in the behavior of the neu- 
tron skin and its bulk and surface contributions along the 
neutron-rich Sn and Pb isotopes can be qualitatively un- 
derstood as follows. Within a major shell, the rms radii 
of the different single-particle orbits are spread around 
their average value. The rms radius of orbits with low 
angular momentum I are larger than the average, while 
the rms radius of orbits with high I are slightly smaller 
but close to the average in the shell. Consequently, the 
outermost region of the density is basically provided by 
the orbits with low I in the last populated major shell. 
On the contrary, orbits with high I in this major shell 
have their most important contribution at shorter dis- 
tances from the center of the nucleus, increasing more 
the bulk than the surface part of the nuclear density. 
Hence, the filling order of the last single-particle orbits is 
crucial to determine if the growing of the neutron radius, 
and consequently the neutron skin, is due to an increase 
of the surface or the bulk of the density. In the case of 
neutron- rich isotopes of Sn above N = 82, the first filled 
orbits are 2/7/2, %P3/2 and 3pi/ 2 . This ordering produces 
an important enhancement of the nuclear surface, which 
can be appreciated from Figs. [7] and El Once midshell 
is filled, the I/19/2 and 1*13/2 orbits start to be apprecia- 
bly populated, increasing the bulk more than the surface 
of the densities and hence of the neutron skin, as can 
be seen in the aforementioned figures. The situation in 
Pb isotopes is just the contrary; above N — 126, the 
first occupied levels are 2g 9 / 2 and lj'15/2, which increase 
the bulk more than the surface. Only near the drip line, 
N = 184, are the low-momentum orbits 3c?5/ 2 , 4si/ 2 , and 
3c?3/ 2 relevant, increasing the surface and quenching the 
bulk contributions to the density and consequently of the 
neutron skin. This behavior in Pb isotopes can also be 
seen in Figs. [7J and |U 

In Figs.[9Uc)-[Hf) we display the results calculated with 
the Skyrme SLy4 force and the NL3 parameter set. The 
bulk and surface contributions qualitatively behave as 
the ones computed with the D1S force discussed previ- 
ously. Of course, some differences are found in comparing 
the results obtained with the different models because of 
their different isovector properties. In models with a stiff 
density dependence of the symmetry energy (for instance, 
NL3), the bulk contribution to the neutron skin is more 
important than in models with a soft symmetry energy, 
as in the case of the SLy4 and D1S forces. This tendency 
can be especially noted in the Pb isotonic chain. 



V. CONCLUSIONS 

Ground-state properties of stable nuclei, such as charge 
radii and binding energies, can be reproduced fairly well 
by using mean-field models with effective interactions 
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such as the Skyrme or Gogny forces and with relativis- 
tic Lagrangians. Although the isoscalar part of these 
models is well constrained, their isovector properties are 
much less determined and the predictions in this sector 
can differ considerably, even for stable nuclei. A typical 
example is the neutron skin thickness of nuclei. For this 
observable theoretical predictions for 208 Pb using non- 
relativistic forces give a value around 0.15 fm while rel- 
ativistic mean-field paramctrizations predict values that 
are almost twice as large. 

The neutron skin thickness of a set of 26 nuclei, dis- 
tributed over the whole periodic table, has been obtained 
from the analysis of experiments with antiprotonic atoms 
[38l |39| combined with the charge radii obtained from 
electron scattering data. One important result of these 
experiments with antiprotonic atoms is that there ex- 
ists a rather clear linear correlation between the neutron 
skin thickness and the overall neutron excess /. The- 
oretical mean-field calculations of the neutron skin also 
show this tendency even more clearly. It is found that the 
relativistic parametrizations systematically predict larger 
neutron skins than the ones computed with nonrelativis- 



tic interactions. This is because the symmetry energy is 
stiffcr in the relativistic models than in the nonrelativis- 
tic models. 

To analyze the experimental data of antiprotonic 
atoms an ansatz of the nuclear densities is needed. The 
two-parameter Fermi distributions have been used often 
to this end. It is found that the experimental data can be 
reproduced by a variety of these Fermi distributions with 
different values of C n — C p and a n — a p . The experimental 
values of the halo factor in 208 Pb are well reproduced by 
distributions with C n ~ C p (halo model) and by Fermi 
densities with both halo and neutron skin (a n rs a p ) con- 
tributions. This latter scenario is also well predicted by 
the nonrelativistic mean-field densities obtained with the 
D1S and SLy4 forces. 

We have also parametrized the mean-field densities 
via two-parameter Fermi distributions. We do this by 
imposing that both mean-field and parametrized densi- 
ties give the same quadratic and quatric moments. This 
parametrization of the mean-field densities also allows 
the neutron skin thickness to be split easily into two con- 
tributions, namely the bulk part and the surface part. It 
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is found that the mean-field neutron skins computed in 
nuclei with I > 0.1 can be shared between non-negligible 
surface and bulk parts. This applies both for stable nu- 
clei investigated in antiprotonic experiments and for drip 
line isotopes in all the theoretical models considered in 
this work. 

To analyze the neutron skin in neutron-rich nuclei, we 
have theoretically studied its variation along the Sn and 
Pb isotopic chains up to the neutron drip line using se- 
lected mean-field models. As expected, AR np shows gen- 
erally linear growing trend with /. However, shell effects, 
which are always present in mean-field calculations, pro- 
duce noticeable departures from this linear dependence 
in nuclei with large neutron excesses. 

Regarding the bulk and surface contributions to the 
neutron skin thickness in Sn isotopes, it can be seen that 
the considered mean-field models point toward more of a 
surface character in stable nuclei. This effect is reinforced 
in the neutron-rich region. In the case of the Pb isotopic 
chain, bulk and surface contributions have similar values 
in stable isotopes, whereas the bulk part is larger than 
the surface part in the more neutron-rich region of this 
isotopic chain. 



where, q — n,p,c denotes the neutron, proton, and 
charge distributions, respectively. The result of Eq. (IA3|) 
can be easily derived by recalling Eqs. dHJ, (fTTj) . (fl"2"T) . and 
([17)1 given in Sec. [TTT1 

From the relation existing between the charge and pro- 
ton rms radii 



0.64 fm 2 



(A4) 



together with normalization condition (|A2[) . one can de- 
rive the parameters a p and C p of the 2pF point proton 
distribution if the parameters a c and C c of the experi- 
mental charge distribution are known (assuming that the 
central density is the same for charge and protons). The 
result for a„ is 



C p / 3Z 



7r y 4irC3p. 



'Or 



and C p is obtained as 



C p — S\ + 02 , 
where Si and S2 are given by the equations 



(A5) 



(A6) 
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Appendix A 

In two-parameter Fermi (2pF) density distributions 

*M = i , ^l° q , (Al) 



+ exp[(r - C q )/a q ] ' 



the number of particles and the mean square radii can 
be approximated, respectively, by the relationships 



N„ 



*C 3 qP0q 1 + q - 



(A2) 



and 



S{ = Tx + J If + T| 



Si. 



n 



with 



Ti = 



21Z 
327rp 0c 



T 2 = - «r 
12 vv 



Ti. 
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(A7) 



(A8) 



(A9) 



Once a p , C p , and the rms radius of the point proton 
density are available, the following relationship between 
th eparameters a n and C n of the neutron distribution can 
be applied: 



r 2\l/2 



>i C] 
7 ir 2 



--f. (A10) 



This expression is obtained by inserting the neutron skin 



thickness Ai?„ 
3| for q = n 



(r 



2\l/2 



IP 



of the nucleus in Eq. 



Therefore, for the same 



„2\l/2 



and AR„ 



values, one has a degenerate family of 2pF neutron den- 
sities depending on one parameter, which can be taken 
to be either C n — C p or a n — a p (recall that C p and a p 

are known). Alternatively, if the values of {r 2 )]/ 2 and 
C n — C p are given, one obtains a family of AR np val- 
ues depending on the parameter a n — a p . As mentioned 
in the main text, the experiments in antiprotonic atoms 
were shown to preferentially support the situation of 2pF 



(A3) density distributions having C n — C p 
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